The Carmeli metric correctly describes 
spiral galaxy rotation curves 
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The metric by Carmeli accurately produces the Tully-Fisher type relation in spiral galaxies, a 
relation showing the fourth power of the rotation speed proportional to the mass of the galaxy. 
And therefore it is claimed that it is also no longer necessary to invoke dark matter to explain the 
anomalous dynamics in the arms of spiral galaxies. An analysis is presented here showing Carmeli's 5 
dimensional space-time-velocity metric can also indeed describe the rotation curves of spiral galaxies 
based on the properties of the metric alone. 
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I. INTRODUCTION 

The rotation curves highlighted by the circular mo- 
tion of stars or more accurately characterised by the 
spectroscopic detection of the motion of neutral hydro- 
gen and other gases in the disk regions of spiral galaxies 
have caused concern for astronomers for many decades. 
Newton's law of gravitation predicts much lower orbital 
speeds than those measured in the disk regions of spi- 
ral galaxies. Most luminous galaxies show slightly de- 
clining rotation curves (orbital speed vs radial position 
from nucleus) in the regions outside the star bearing disk, 
coming down from a broad maximum in the disk. Inter- 
mediate mass galaxies have mostly nearly flat rotation 
speeds along the disk radius. Lower luminosity galaxies 
usually have monotonically increasing orbital velocities 
across the disk. (See for a good review.) The tradi- 
tional solution has been to invoke halo 'dark matter' p] 
that surrounds the galaxy but is transparent to all forms 
of electromagnetic radiation. In fact, astronomers have 
traditionally resorted to 'dark matter' whenever known 
laws of physics were unable to explain the observed dy- 
namics. 

In 1983 Milgrom introduced his MOND [HEEl, 
an empirical approach, which attempts to modify New- 
tonian dynamics in the region of very low acceleration. 
Newton's law describes a force proportional to r -2 , where 
r is the radial position,but Milgrom finds that ar" 1 law 
fits the data very well [lj . Others have also attempted to 
formulate modified force laws, such as Disney Q, Wright 
[r?| and Carmeli [HEl- The latter formulated a modifica- 
tion and an extension of Einstein's general theory, in an 
expanding universe taking into account the Hubble ex- 
pansion, which imposes an additional constraint on the 
dynamics of particles 0- 

Carmeli believes the usual assumptions in deriving 
Newton's gravitational force law from general relativity 
are insufficient, that gases and stars in the arms of spiral 



galaxies are not immune from Hubble flow. As a con- 
sequence a universal constant ao (in this case, slightly 
different to Milgrom's) is introduced as the minimum ac- 
celeration in the cosmos. 

Using this theory Carmeli Q successfully provided a 
theoretical description of the Tully-Fisher law, deriv- 
ing the proportionality between the fourth power of the 
galaxy's rotation speed (v) and its mass (M), vis-a-vis 



v 4 = -GMa , (1) 



where G is the Newtonian gravitational constant. 
Equation JIJ can be re-written as v 4 = (^r-frCio) 

and therefore the positive square root is v 2 jr - 
where is the Newtonian gravitational acceleration. 
Hence the latter is consistent with Milgrom's phenomeno- 
logical approach in the low acceleration limit. 

II. CARMELI'S 5D COSMOLOGY 

In the weak gravitational limit, where Newtonian grav- 
itation applies, it is sufficient to assume the Carmeli met- 
ric with non-zero elements <?oo = l + 0/c 2 j 544 = l+'vV' 7 " 2 ; 
gkk = — 1; [k = 1,2,3) in the lowest approximations 
in both 1/c and 1/t. Here a new constant, called the 
Hubble-Carmeli constant, is introduced r « 1/Hq . The 
potential functions <f> and ip are determined by Einstein's 
field equations and from their respective Poisson equa- 
tion, either 
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V 2 </> = 4nGp m or V 2 ?/> = 



AttG 



(2) 
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where p m is the mass density and ao the universal char- 
acteristic acceleration ao = c/r. As usual c is the speed 
of light in vacuo. 

The Hubble law describes the expansion of the cos- 
mos and the matter embedded in it. Therefore the 
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line element for any two points in this 'new' space-time- 
velocity is ds 2 — gooC 2 dt 2 + g^dx k + g^r 2 dv 2 = 0. 
The relative separation in 3 spatial coordinates r 2 = 
(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 and the relative velocity between 
points connected by ds is v. The Hubble-Carmeli con- 
stant, r, is a constant for all observers at the same epoch, 
therefore may be regarded as a constant on the scale of 
any measurements. 

The equations of motion (B.62a and B.63a from HJ) to 
lowest approximation in 1 / c are reproduced here, 



d 2 x k 



1 8<f> 
2dx^' 



(3) 



This is the usual looking geodesic equation derived from 
general relativity but now in 5 dimensions. And the sec- 
ond is a new phase space equation derived from the the- 
ory 



r is the familiar Newtonian equation. The first term on 
the right-hand side of J7|) Carmeli neglects as small on a 
galaxy scale. But that would only be true if the matter 
density (Q) is a descriptor of space curvature on a much 
larger scale than the galaxy. If it describes the local 
density then the term is not insignificant. 



B. Post-Newtonian circular motion 
From (0J and l|5b"|l it follows that 



d 2 x k 
dv 2 
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r 



Integrating with respect to dx k © becomes 



(8) 
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dv 2 
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2dx k ' 



(4) 



The solution of Einstein's field equations in 5D result in 



o 2GM 
T^r 2 , 



i - n 



2GM t 2 



(5a) 
(5b) 



where fl is the matter density expressed as a fraction of 
the critical or 'closure' density, which in this model is 
defined by p c = 3/8ttGt 2 . 
From (0 and (J^U it follows 



d 2 x k 



,,,, - l -^(r% k+ GM(l r ),. 



A. Post-Newtonian force law 

By carrying out the differentiation with respect to r, 
© becomes a new post-Newtonian force equation 



g( r ) = - 



i - n 



GM 



(7) 
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GMt 2 1 



(9) 



which is the new post-Newtonian equation for circular 
motion. 

For r small, © becomes equation (8) of Because 
the new dimension (u) in the Carmeli metric is con- 
structed as an analogue to the time co-ordinate in the 
usual 4D spacetime, the meaning of the new equation is 
identified by the substitutions v — > t, t — > c. In this 
case, with these substitutions, and Q = 1 representing 
Euclidean space, we recover the usual equation for circu- 
lar motion. 



C. Ontology 

A possible ontology for @ is that even though galax- 
ies are constrained by gravity against expansion in the 
radial direction, they are free to expand tangentially to 
the radius i.e. azimuthally in the plane of the galaxy. 

It is known that some galaxies have portions of their 
disks, and the extended gas regions bey ond, that rotate 
in opposite directions to each other [l5j. From the sug- 
gested ontology, it would be expected that this situation 
could occur. As a new galaxy develops the direction of 
rotation is determined by the initial angular momentum 
of the system. For proto-galaxies with low initial angu- 
lar momentum different portions of the developing disk 
can rotate in response to the Hubble flow force in the az- 
imuthal direction, some parts prograde, some parts ret- 
rograde. 



For f2 less than critical (fi < 1) in or near a galaxy this 
means an additional force inwards is applied to the test 
particle. For fl more than critical (0 > I) it represents 
an additional outward force. The solution of Q for small 



D. Phase space equation 

After integrating (0 and solving for radial distance 
(r) as a function of velocity (u) with the condition v(r = 



3 



0) = 0, we get 



For Q >• 1, © becomes 



/2GMr 2 

V n- i 



sinh 3 




(10) 



which, to lowest order in v/c, reduces to 



/GMt 2 



(11) 



This is equation (9) of Q . For typical-mass galaxies l(TT|) 
describes orbital speeds of hundreds of kilometres per 
second. On small r scales, or scales of small local orbital 
motion, the effect of is small. Assuming r small 
so the first term is negligible and integrating (JJJ) with 
respect to dr we get the familiar equation for circular 
motion 



GM 



(12) 



where v is the circular or tangential velocity in the sys- 
tem. Simultaneously solving (|T^|l and (fTO by eliminating 
r yields a Tully-Fisher type relation as in J3>. This is the 
result of Q. 



III. GALAXY ROTATION CURVES 



The first term in (|5a|l and hence J7|) is valid where 



Ic 1. By substituting p c into 



and assuming 



51 ^ 1 it can be shown this condition is satisfied when 
f 2 Pm < 10 27 (in SI units), where p m is the local bary- 
onic matter density. For radial distances within a galaxy 
(r) on the scale of a few kpcs (1 kpc « 3 x 10 19 m) the 
density of material p m < 1CP 14 g cm~ 3 . This is true for 
most regions in a galaxy; the exception would be near a 
compact object. 

By rearranging (J7J) using £1 = p m j p c where ft ^ 1 we 
get 



AnGpr, 
3a 2 , 



GM 1 



(14) 



Notice ao (= c/t) explicitly appears in this equation. In 
the limit of zero mass and zero distance (0 can be re- 
written as 



s) 



(15) 



which is the Hubble relation. Each term on the rhs of 
(|14fl contributes half the gravitational radius cancelling 
the radius in the denominator. 

Now (ITU) , like (fH^ . is galaxy specific and depends on 
the matter density in the disk. Using (|13f) a new expres- 
sion for the speed (v) of circular motion of stars and gases 
will be derived as a function of distance (r) from the cen- 
tre of a galaxy. Real values of typical densities will be 
used to see if a connection can be found to the Milgrom 
formula, which means g(r) — > ^/g^ao as r — - > oo. In this 
case it turns out ao is different by a constant factor. 

By substitution of O — 1 — > 8ttGt 2 p m /3 into (|l()(l . for 
fl » 1, we get 



/ 3M 



\A-Kpr, 



sinhi QV^GWni- 2 ) . (16) 



By integrating l|13|) with respect to dr we obtain the new 
equation for circular motion 



4nGp 



m 2 

r 



GM 



(17) 



for £1 S> 1. Notice this equation contains a new term, 
which adds additional velocity to the stars and gases cir- 
culating in the disk region of the galaxy. But for small r 
the term is insignificant. 

If we initially assume that the argument of the hyper- 
bolic sine function in l|ltj|) is much less than unity, that 
is, 



9(r) 



8nGp 71 



GM 



(13) 



Note there is no term dependent on the Hubble-Carmeli 
time constant (r). Equation (|13fl is a total local force 
law, and a post-Newtonian equation. It needs significant 
investigation on the scale of galaxies, clusters etc. It 
indicates a gravity shielding type effect by the matter 
density field. For small radii from a central gravitating 
body, the first term is insignificant with respect to the 
second term and (|13fl becomes the familiar Newtonian 
equation. For the Sun-Earth system the second term is 
about 30 orders of magnitude larger than the first. 



-\/3nGp m T 2 =i«l 
c 



(18) 



= 10 23 g cm 3 . 



which implies v <C 150 km s _1 for p ri 
This is reasonable in some disk regions of galaxies, but 
when higher velocities are measured this approximation 
cannot be used. 

Using the expansion sinhi (x) « x 5 + x s /9 for 
collecting terms lowest in x and substituting r from ll(j|) 
into (|17fl for v <^ c for circular motion, we get 



= \ -GMa 1 



3irGp 71 



(19) 
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For very low matter density (p m — > 0) this becomes 
the Tully-Fisher type law for galaxies as in (JIJ. However 
where the approximation ljl8jl is true, it also follows that 

3 -^v^(lV^G^) 2 <<l. (20) 
Hence l|19[) becomes 

and the solution for circular velocity is simply found. 

« 2 =\/I^( i - ? ^HP^) 1 > 22 > 

Equation (|22[) applies in all regions of spiral galaxies 
with densities that 10~ 30 g cm' 3 < p m < 10' 23 g cm' 3 . 
Densities in this range yield circular velocities in the 
disk regions like 85 km s' 1 for galaxies with mass 
M = 10 9 M Q and 150 km s _1 for galaxies with mass M = 
10 10 M Q (where M Q is a solar mass unit — 2 x 10 30 kg). 
For larger masses, the higher velocities mean the approx- 
imation is invalid. Therefore for higher velocities and/or 
densities p m > 10' 23 g cm' 3 ((TCI) must be used without 
approximation. 



V. TOY MODEL: DENSITY PROFILE 

Let us assume a simple density distribution for a spiral 
galaxy valid in the disk region, expressed in cylindrical 
co-ordinates (r, 0, z) but with no (^-dependence on mass 
density 

p m (r,z) = Ar- a e- 0z , (26) 

where A, a and (3 are parameters to be determined. Here 
(3 = 1/zq, and zq is an exponential scale factor of order 
of several pc for a galaxy. Therefore in this simple toy 
model the accumulated mass of the galaxy M(r) can be 
written as the sum of two parts, the bulge mass (Mo) 
plus the disk mass, hence 

M(r) = M + 2vr J Q r rdr f* z dz p rn (r, z) 

= M + 2ttA J Q r r^dr 2 /„* e' z ' z °dz, (27) 

which yields mass as a function of r. 

M(r)=M + ^V-«, (28) 
2 — a 

where q^2. 

Also assume the density is uniform in z and valid for 
z = 0, hence in the appropriate units, 

Pm (r) = Ar' a (29) 



IV. MILGROM'S MOND 

We can write gAi(r) — v 2 /r where Qm{t) is the 
Milgrom acceleration inferred from circular motion and 
therefore it follows from 1221) that 

9m(t) = \j^gN{r)a ^1 - ^a"' \j\ GM ( r ) a ^j - 

(23) 

where M(r) is the enclosed mass at radius r. By squaring 
it follows 

9Ur) = 9N(r)a' (l - H^^ G M(r)a'^ , (24) 

where gAr(r) is the Newtonian acceleration and a' = 
|ao ~ 5.2 x 10 -10 m s~ 2 , depending on the precise value 

Of T. 

In the region of low acceleration on the edge of a galaxy 
where p m — + as r — * large, and M{r) — > a constant, l|24() 
becomes 

9M( r ) = 9N(r)a' (25) 

which is the low acceleration limit of Milgrom's phe- 
nomenological law. 



Substituting l(77|) and into J22J) we get for circular 
motion 

Assuming zq — lpc, (|30|l was fit to the velocity- 
distance data of NGC6503 from [jj] with the results shown 
in curve 1 (dashed line) in figure ^ The best fit param- 
eters for Mo, A and a were determined. The fit yielded 
a first order density dependence on r. It also showed 
a good agreement, showing the correct trend, especially 
in the disk region, considering the simple model used. 
The total mass of the galaxy is then M(r = 15 kpc) = 
M + 4nAz r = (1.836 + 0.165) x 10 9 M . Therefore 
the disk comprises 8.9% of the total galaxy mass and 
the total mass of the galaxy equals 2.0 x 10 9 M©. This 
compares with published value of 4.8 x 10 9 M Q for 
the luminosity mass of the galaxy. It is underestimated, 
which is the reverse of the dark matter problem. 

But most importantly no dark matter is needed to ex- 
plain the rotation curve. The shape of the derived de- 
pendence of the tangential speed of the stars to their 
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FIG. 1: Fit to the velocity-distance data of the spiral galaxy 
NGC6503. Filled circles are the data, and the dashed line 
(curve 1) the fit from 13111 . The parameters Mo = 1.836 x 
1O 9 M and a = 1.0075 were determined from the fit. Curve 
2 (solid line) represents the fit using an undulating density 
function as shown in figure |3] 



radial distance from the centre of the galaxy is in accor- 
dance with expectation. A more comprehensive model 
though is needed especially for the region close to the 
central bulge. That should involve a gas component as 
well as a more precise model of the enclosed mass. Also 
the approximate model used here would not apply as the 
density and orbital speeds will be greater in many cases. 
However the purpose of the toy model presented here is 
to show that it can work in the outer disk region of a 
galaxy and realistic though much smaller masses result 
from the fit than when halo dark matter is invoked. 



VI. GENERAL SOLUTION 

After substituting l|16(l into l|17fl . eliminating r, a tran- 
scendental equation of the following form results, 
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FIG. 2: Solution of I31H with orbital velocity (v) on the 
vertical axis in km s _1 , the galaxy mass (M) in 10 9 Mq units 
and the matter density (p m ) in g cm~ 3 
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(31) 



FIG. 3: Simulated density fluctuations (in g cm 3 ) as a 
function of disk radius (in kpc) 



where 8 = (3M/47rp m ) 1/3 , ? = y/3nGp m /a , 7 = GM. 

Equation (|31|l must be solved for the motions of stars 
and gases in any region of a galaxy. An exact solution 
of |^nj is of course not analytically possible but the so- 
lution may be visualized using a 3D plot generated by 
the Mathematics software package. This is shown in fig- 
ure It is apparent that as we follow a line of constant 
mass (M) we see the tangential velocity (v) increase as a 
function of local matter density (p m )- The more massive 



the galaxy the less sensitive the curves are to changes in 
density. 



A. Density fluctuations 

As is noted in many galaxy rotation curves the velocity 
often rises and falls as a function of radial distance. This 
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TABLE I: Expected mass of Milky Way galaxy 
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FIG. 4: Various rotation curves generated from l.'ilt . See the 
text for details 



a fig.4 ofJTal 

see ref lltl . Note, there is so me evidence for a Keplerian declining 
rotation curve beyo nd 1 7 kpc jlflfl . The mass distribution > 22 kpc 
is still controversial Il5l. 



galaxy with a density dependence of r 8 (i.e. a = 0.8) 
and M = 2.8 x 10 9 Af© but with 14.4 x 10 9 A/ Q in the 
disk out to 15 kpc. That is 83% in the disk. This curve 
shows a continuous rise as a function of distance as most 
of the mass is exterior to the nucleus. 



B. Galaxy mass 



is indicative of changes in density in the disk region. In 
order to simulate the undulating density fluctuations in 
the arms of a spiral galaxy, which commonly are of the 
order of a few tens of kilometres per second , I have 
chosen a density function as shown in figure [3] This 
function has a mean r _1 dependence on radius as per the 
result of the toy model above. Using this function and 
a mass derived from l|28l) , replacing the density function, 
(|31|l was solved for 20 fixed points along the radius, where 
mass and density are both functions of radius. The mass 
Mq was adjusted to get the best fit to the simulated data 
of NGC6503 and the result is shown in curve 2 (solid line) 
in figure ^ The mass of the galaxy was then determined 
to be 3.2 x 1O 9 M with 14% of the mass in the disk out 
to 15 kpc. This is somewhat different to the toy model 
result above, and the fit is not as good, however I have 
tried to simulate the effect of density fluctuations. Note 
it has not been attempted to simulate density and mass 
within the central 5 kpc. That will be left to future work. 

By changing the mass of the central nuclear bulge from 
Mo = 2.8 x 10 9 M Q by 1 and 2 orders of magnitude the 
following curves in figure^Jwere generated. Curve 1 is the 
same curve as curve 2 in figure ^ curve 2 is with Mq = 
28 x 10 9 M Q and curve 3 is with M = 280 x 10 9 M Q in 
the central core. Notice the stronger density fluctuations 
and the rigid body steep nuclear rise. Normally the latter 
would fall well within a kpc but I have not attempted to 
model the nuclear bulge region. 

In each of the above the density dependence was r^ 1 
(i.e. a = 1) and with about 10% of the mass in the disk 
out to 15 kpc, as in the toy model above. By altering the 
density fall off we get a rising velocity curve. Curve 4 of 
figure 01 represents a galaxy with a 1O 1O M core and only 
7.5% in the disk out to 15 kpc and a density function 
falling as r~°- 9 (i.e. a = 0.9). Curve 5 represents a 



Equation l|31|) can be rearranged to calculate the mass 
of a galaxy from the known velocity and an estimate on 
the density. 



M 



1/2 T • i/ \ 

v sinn(c,v) 
cosh 3 (sv) 



(32) 



Equation (|32(l is a post-Newtonian equation and con- 
tains, besides constants, a measured velocity and the 
matter density which needs to be estimated in the region 
where the velocity is measured. The universal accelera- 
tion constant clq is contained in <;. Here it was assumed 
that a = 7.8 x 10~ 10 m s~ 2 after Carmeli 0. 

Using the measured orbital speeds around the Milky 
Way galactic centre we can calculate an enclosed mass 
from H32[) and compare it to the enclosed mass derived 
from the Newtonian circular motion equation l|12fl as well 
as the commonly cited mass of 10 n M Q for the Galaxy. 
Table [I] shows the results of the calculations. Our sun's 
position is within the 9.95 kpc and the enclosed mass is 
4 times smaller from the post-Newtonian equation when 
a reasonable density is assumed. If a density 2 orders of 
magnitude larger is assumed, the mass is 7 times smaller. 
At 18.4 kpc the enclosed mass is 3 times smaller depend- 
ing on the matter density. 

From Table |U is can be clearly seen that the 'missing' 
mass or the dark matter that is supposed to haunt the 
halo regions of our galaxy is the result of the incorrect 
equation of circular motion. Depending on the matter 
density distribution through the galaxy the total mass 
may in fact be at least 4 times less than previously cal- 
culated. 
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VII. CONCLUSION 

Carmeli successfully predicted the accelerating uni- 
verse 3] two years before the announcements 8,14]. His 
new metric has validity on the scale of the universe with- 
out assuming any dark matter 0- Here it is also shown 
that it may be the solution to the rotation curve anomaly 
in the outer regions of spiral galaxies. 

Equations of motion have been derived from Carmeli's 
metric, which produce a Tully-Fisher type relation and 
describe the rotation curves in spiral galaxies without the 



need for non-baryonic halo dark matter. New equations 
for circular motion are discovered. A theoretical compar- 
ison is made with Milgrom's MOND phenomenology and 
some agreement found in the low acceleration limit. 

Based on the 5D space-time-velocity of the Carmeli 
metric the assumed dark matter in galaxies to explain 
the anomalous rotation curves is no longer needed. This 
situation is similar to the case when Einstein introduced 
his general theory the advance of the perihelion of Mer- 
cury was explained adequately without the need for dark 
matter as had earlier been believed. 
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